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Abstract 

We obtain a discrete time analog of E. Noether's theorem in Opti- 
mal Control, asserting that integrals of motion associated to the discrete 
time Pontryagin Maximum Principle can be computed from the quasi- 
invariance properties of the discrete time Lagrangian and discrete time 
control system. As corollaries, results for first-order and higher-order dis- 
crete problems of the calculus of variations are obtained. 
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1 Introduction 

Most physical systems encountered in nature exhibit symmetries: there exists 
appropriate infinitesimal-parameter family of transformations which keep the 
system invariant. From the well-known theorem of Emmy Noether |261 127| . 
one can discover the integrals of motion from those invariance transformations. 
Noether's theorem plays a fundamental role in modern physics, and is usually 
formulated in the context of the calculus of variations: from the invariance 
properties of the variational integrals, the integrals of motion of the respective 
Euler-Lagrangc differential equations, that is, expressions which are preserved 
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along the extremals, are obtained. The result is, however, much more than a 
theorem. It is an universal principle, which can be formalized in a precise state- 
ment, as a theorem, on very different contexts and, for each such context, under 
very different assumptions. Let us consider, for example, classical mechanics or, 
more generally, the calculus of variations. Typically, Noether transformations 
are considered to be point-transformations (they are considered to be functions 
of coordinates and time), but one can consider more general transformations 
depending also on velocities and higher derivatives |1U| or within the broader 
context of dynamical symmetries (201 ■ For an example of an integral of mo- 
tion which comes from an invariance transformation depending on velocities, 
see [23|. In most formulations of Noether's principle, the Noether transforma- 
tions keep the integral functional invariant (cf. e.g. §1.5]). It is possible, 
however, to consider transformations of the problem up to an exact differen- 
tial (cf. e.g. [SU p. 73]), called a gauge- term [j^l • Once strictly- invariance 
of the integral functional is no more imposed, one can think considering addi- 
tional terms in the variation of the Lagrangian - see the quasi-invariance and 
semi-invariance notions introduced by the author respectively in |38j and |39j . 
Formulations of Noether's principle are possible for problems of the calculus of 
variations: on Euclidean spaces (cf. e.g. [15)1 or on manifolds (cf. e.g. dJ); 
with single or multiple integrals (cf. e.g. 5 ); with higher-order derivatives (cf. 
e -g- D); with holonomic or nonholonomic constraints (cf. e.g. 021 Ch. 7], 
|35|1; and so on. Other contexts for which Noether's theorems are available 
include supermechanics [SJ, control systems 02 EH an d optimal control (see 
e.g. El E3 EH!)- For a survey see [123 Here we are interested in pro- 
viding a formulation of the Noether's principle in the discrete time setting. For 
a description of discrete time mechanics, discrete time calculus of variations, 
and discrete optimal control see, e.g., ^3 El EE1 Ell , HI) an( A ^21- Illustrative 
examples of real-life problems which can be modeled in such framework can be 
found in ^21 Ch. 8]. Versions of the Noether's principle for the discrete calculus 
of variations, and applicable to discrete analogues of classical mechanics, ap- 
peared earlier in |161 1171 1211 [31 ^2 ^] , motivated by the advances of numerical 
and computational methods. There, the discrete analog of Noether's theorem 
is obtained from the discrete analog of the Euler-Lagrange equations. To the 
best of our knowledge, no Noether type theorem is available for the discrete 
time optimal control setting. One such formulation is our concern here. The 
result is obtained from the discrete time version of the Pontryagin maximum 
principle. As corollaries, we obtain generalizations of the previous results for 
first-order and higher-order discrete problems of the calculus of variations which 
are quasi- invariant and not necessarily invariant. 



2 Discrete-Time Optimal Control 

Without loss of generality (cf. [^ §2]), we consider the discrete optimal control 
problem in Lagrange form. The time k is a discrete variable: k s Z. The horizon 
consists of N periods, k = M, M+l, . . . ,M + N — 1, where M and N are fixed 
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integers, instead of a continuous interval. We look for a finite control sequence 
u(k) e W, k = M,...,M + N — 1, and the corresponding state sequence 
x(k) £ K™, k = M, . . . , M + N, which minimizes or maximizes the sum 

M+N-l 

J[x(-),u(-)] = E L(k,x(k),u(k)) , 

k=M 

subject to the discrete time control system 

x{k + l) = ip{k,x{k),u{k)) , k = M,...,M + N-l, (1) 

the boundary conditions 

x{M) = x M , x{M + N) = x M +N , (2) 

and the control constraint 

u(k) 6S]Cl r , k = M, . . . , M + N - 1 . 

A sequence-pair (x(k), u(k)), k = M, . . . , M + N — 1, satisfying the recurrence 
relation JQ) and conditions |[5J, is said to be admissible: x(k) is an admissible 
state sequence and u(k) an admissible control sequence. Functions L(k, x, u) : 
{M, ... ,M + N —1} x K" x M r — » M and ^(fc, x, u) : {M, . . . , M + N - 1} X 
M™ x R' r — > R n are assumed to be continuously differentiable with respect to a; 
and u for all fixed k = M, . . . ,M + N — 1 , and convex in u for any fixed k and 
x. They arc in general nonlinear. The control constraint set f2 is assumed to be 
convex. The problem is denoted by (P). 

Remark 2.1. For continuous optimal control problems, the convexity assump- 
tions we are imposing are not needed in order to derive the Pontryagin maximum 
principle |3(Jj . This differs from the discrete time optimal control setting. Our 
hypothesis can be, however, weakened to directional convexity or even more 
weak conditions (see [21], [121 §8-3] and references in |321 Ch. 6] and [24)1 . 

Remark 2.2. It is possible to formulate problem (P) with the first-order differ- 
ence equations i n terms of the forward or backward difference operators A 
or V, defined by Ax(fc) = x(k + 1) — x(k), Vx(k) = x(k) — x(k — 1). The results 
of the paper are written in those terms in a straightforward way. 

The following theorem provide a first-order necessary optimality condition 
(cf. e.g. 0J §3.3.3], [22], |321 Ch. 6]) in the form of Pontryagin's maximum 
principle |3(J|. For a good survey on the history of the development of maximum 
principle to the optimization of discrete time systems, we refer the reader to |24| . 

Theorem 2.1 (Discrete-Time Maximum Principle). If (x(k),u(k)) is 
a minimizer or a maximizer of the problem (P), then there exists a nonzero 
sequence-pair (i/'o, ip(k)) 7 k = M + 1, . . . , M + N , where ipo is a constant less or 
equal than zero and ip(k) G M™, such that the sequence-quadruple 

(x(k),u(k),i[> ,i>(k + l)) , k = M,...,M + N-l, 

satisfies: 
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(i) the Hamiltonian system 

x(k + l) = ^(k,x(k),u(k),%,^(k + l)) , k = M,...,M + N-l, 
i/>(k) = §f (k, x(k),u(k), ifo, + 1)) , k = M + 1, . . . , M + N - 1 ; 

(3) 

(ii) the maximality condition 

H (k, x(k),u(k), ip , ip(k + 1)) = maxH (k, x(k),u, ipo, ip{k + 1)) , (4) 

k = M, ...,M + N -1; 
with the Hamiltonian 

H (k, x, u, ipo, ip) — ipoL(k, x,u) + if> ■ (p(k, x, u) . 

Remark 2.3. The first equation in the Hamiltonian system is just the control 
system Q). The second equation in the Hamiltonian system is known as the 
adjoint system. The multipliers ip{-) are called adjoint multipliers or co-state 
variables. 

Remark 2.4. In the absence of the initial conditions x(M) = Xu and/or terminal 
conditions x{M + N) = xm+n, there corresponds additional conditions in the 
Discrete-Time Maximum Principle called transversality conditions. Our version 
of Noether's theorem only require the use of the adjoint system and maximality 
condition. Therefore, the result is valid under all types of boundary conditions 
under consideration. 

Definition 2.1. A sequence-quadruple (x(k), u(k), t/'oj ip(k + 1)), k = M, . . . , M+ 
N — 1, ipo < 0, satisfying the Hamiltonian system and the maximality condi- 
tion, is called an extremal for problem (P). An extremal is said to be normal if 
ipo 5^ and abnormal if ipo = 0. 

Remark 2.5. As we will see on Section|5| there are no abnormal extremals both 
for first-order and higher-order discrete problems of the calculus of variations. 
In particular, there are no abnormal extremals for problems of discrete time 
mechanics. For our general problem (P), however, abnormal extremals do exist. 
In fact, they happen to occur frequently. For a throughout study of abnormal 
extremals see 

3 Integrals of Motion 

We obtain a systematic procedure to establish integrals of motion, i.e., to es- 
tablish expressions which are preserved on the extremals of the discrete optimal 
control problem (P), from the (quasi-) invar iance properties of the discrete La- 
grangian L {k, x(k), u{k)) and discrete control system x(k+X) = if (k, x(k), u(k)). 
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Definition 3.1. Let X : {M, . . . ,M + N — 1} x R n x £1 x S(0;e) -> R", 
£ > 0, B(0;s) = |s = (si, . . . ,s p )| ||s|| = v / J2i=i( s i) 2 < £ }> be an infinitesi- 
mal p-parameter transformation such that for each k, k = M, . . . , M + N — 1, 
X(fc, •, •, •) is continuously differentiable with respect to all arguments, and such 
that X(k, x, u, 0) = x for all k = M, . . . , M + N - 1, x G R", and it G O. If 
there exists a real function $ (fc,x,it, s) and for all s G £>(0;e) and admissible 
(x(k), u(k)) there exists a control sequence u(fc, s), 0) = u(fc), such that: 

L (k, x(k), u(k)) + A$ (k, x(k),u(k), s) + 6 (k, x(k), u(k),s) 

= L (k, X (k, x(k),u(k), s) ,u(k, s)) , (5) 



X(k + l, x(k + 1), u(k + 1), s) + 6 0, x(k),u(k), s) 

— ip (k, X (k,x(k),u(k), s) ,u(k, s)) , (6) 

for each k = M, . . . , M + N — 1 and where S(k, x, u, s) is an arbitrary function 
satisfying 

dS(k,x,u,s) 

=0, i=l,...,p, (7) 

° S i s=0 

for each k, x, u, then the problem (P) is said to be quasi-invariant with 
respect to the transformation X(Jk, x, u, s) up to the difference gauge term 
$ (k,x(k),u(k), s). 

Remark 3.1. In the relation (JSJ), A is the forward difference operator: 

A<E> (k, x(k),u(k), s) = $ (k + 1, x(k + 1), u(k + 1), s) - $ (fc, x(k), u(k),s) . 

Remark 3.2. When 6 = and $ = 0, we have (strict- )invariance. The term 
quasi-invariant refers to the possibility of S to be different from zero. 

Theorem 3.1 (Discrete-Time Noether Theorem). If(P) is quasi-invariant 
with respect to the p-parameter transformation X up to the difference gauge term 
in the sense of Definition UTll then all its extremals (x(k), u(k), ip(k)), 
k = M, . . . , M + N — 1, satisfy the following p expressions (i = 1, . . . , p): 



d 

osi 



+ 4>{k)- ^-X{k,x(k),u(k),s) 



= constant. 



Remark 3.3. The integrals of motion obtained by Theorem 13. II are "momen- 
tum" integrals. Due to the fact that time k is discrete, one can not vary k 
continuously and, for that reason, one can not obtain the "energy" integrals as 
in the continuous optimal control case (cf. (SZHHIj)- To address the problem 
another method needs to be developed. This will be addressed in a forthcoming 
paper. 

Remark 3.4. Together with the continuous results in |371 I4l)j. Theorem 13.11 
provides a framework to obtain a generalization of Noether's theorem for hybrid- 
systems. This and related questions are under study and will be addressed 
elsewhere. 
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Proof. Let (x(k),u(k), tpo,ip(k)) be an extremal for problem (P). Differentiating 
J5j and Jfjjl with respect to the parameter s^, i = l,...,p, and setting s 



(si, . . 
«(*)): 



0, we get (recall J7J) and that X (&;, u(k), 0) = a;(fc), u(fc, 0) 



A — $(fc,3j(k),u(fc),«) 

ds,; 



9i d 

= — (k,x(k),u(k))- —X(k,x(k),u(k),s) 

dL d 
+ -g^ (k, x(k), u(k)) ■ —u (k, s) 



s=0 

(8) 



X (k + l,x{k + l),u(k + 1), s) 



s=0 



^ (fc, x(k), u(k)) ■ (fc, x(k), u(k),s) 



s=0 



+ ^ (k, x(k), u(k)) ■ -^-u (k, s) 



(9) 



From the adjoint system tp(k) = ^ (k, x(k),u(k), ipo, ip(k + 1)), we know that 

-Vo|^ (k, x(k),u(k)) = ^(fc + 1) • || (fc, u(fc)) - V(*0 , 
and multiplying (JSJ) by — one obtains: 



% ( A —§(k,x(k),u(k),a) 



dL d 
- — (k, x(k), u(k)) ■ —u (k, s) 
ou osi 

- (<•(/•' \ )-^(k,x(k),u(k))-ij(k)) ■ JLx(k lX (k),u(k),s) 



= 0. 



(10) 



As far as u(k, 0) = u(k), according to the maximality condition of the Discrete- 
Time Maximum Principle the function 

s i > ipoL (k, x{k), u{k, s)) + ip(k + 1) ■ <p (k, x(k), u{k, s)) 

attains its maximum for s = 0. Therefore, 

8 



ds, 



(t^oL (k, x(k),u(k, s)) + ip(k + 1) • <p (k, x(k),u(k, s))) 



= 0, 



that is, 



dL d 

i(>Q—(k,x(k),u(k))- —u(k,s) 
ou osi 



+ ip(k + 1) ~ (k, x(k), u(k)) ■ -^u(k, s) 



0. (11) 



G 



From l(TU|l and ifTTJl it comes 



d 

i> A —<I>(k,x{k),u(k),s) 



dip 



s=0 



f ( iP(k + l)--^(k,x(k),u(k))-i>{k) 



ds 



■X (fc, x(k), s) 



■%l>{k + 1) ~ (fc, x(k),u(k)) ■ ^-u(k, s) 

OU OS; 



s=0 
= 0. 



s=0 



Using jnj, this last equality is equivalent to 



A (V>o -^-${k,x(k),u(k),s) 
The proof is complete. 



s=0 



s=0 



= 0. 



□ 



4 An Example 

We now illustrate the use of Theorem 13.11 by the following example (n = 3, 
r = 2, ft = M 2 ): 

^ (ui(fc)) 2 - (u 2 (k)) 2 — > extr , 

fc 

Xi(k + 1) = x 2 (k) + ui(k) , 
x 2 {k + 1) = x\(k) + u 2 (k) , 
^xa(k + 1) = x 2 (k)u\(k) , 

subject to fixed endpoints. In this case the Hamiltonian is given by 



H(xi,X 2 ,Ui,U 2 , 1p , 1pl,1p 2 , 4>3) 

= ip L(ui,u 2 ) +ip 1 ip 1 (x 2 ,u 1 ) +?p 2 <p 2 (xi,u 2 ) + ip 3 ip 3 (x 2 ,ui) , 

with L(ui,u 2 ) = u\ - w|, ipi{x 2l u\) — x 2 + ui, f 2 (xi,u 2 ) = xi + u 2 , and 
ip 3 (x 2 ,ui) = x 2 u\. From the adjoint system we get the evolution equations 

Mk) = Mk + 1) + Mk + l)ui(k) , 
Mk) = 0, 

while from the maximality conditions we get (i^^ = 0) 

Mk + l) = -2ipoui(k), 
Mk + I) = 2^u 2 {k) . 
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There are no abnormal extremals for the problem, and one can fix = — h- 
The extremals are obtained solving five difference-equations of order one, 

' xx{k + 1) = x 2 (k) + ipx{k + 1) , 

x 2 {k + l) = xi(fc) - ip 2 (k + l) , 
< x 3 (k + l) = x 2 (k)t/j 1 (k + 1) , 

^(fc + l) = ^ 2 (fc), 
> 2 (fc + l) =Mk), 

together with the boundary conditions (or the transversality conditions), by 
standard techniques. On the other hand, the problem is quasi-invariant with 
respect to the one-parameter (p = 1) transformations 

Xi (xi{k),s) =x-i(k) + 2s, 
X 2 (x 2 (k),s) = x 2 (k) + s , 
X 3 (xi(k),x 3 (k),s) = x 3 (k) + sxi(k) , 

up to the difference gauge term <!> (cci(fc), x 2 (k), s) — 2 (xi(k) + x 2 (k)) s. To see 
that we choose 

Ui(k, s) = Ui(k) + s , u 2 (k, s) — u 2 (k) - s , 

in the Definition 13. f I We notice that X\ (xi(k),0) = xi(k), X 2 (x 2 (k),0) = 
x 2 (k), X 3 (xi(k),x 3 (k),0) = x 3 (k), u 1 (k,0) = Ui(k), and u 2 (k,0) = u 2 {k). 
Direct verifications show that the quasi-invariance conditions are satisfied: 

L( Ul (k, s),u 2 (k, a)) = (ui(k)) 2 - (u 2 {k) f + 2 ( Ul (k) + u 2 {k)) s 

= L (u 1 (k),u 2 (k)) + 2 ( Xl (k + 1) - x 2 (k) +x 2 (k + 1) - x^k)) s 
= L( Ul (k),u 2 (k)) + A$(a;i(fc),a; 2 (fc),s) , 

<Pi (X 2 (x 2 (k), s) , ui(k, s)) = x 2 (k) + ui(k) +2s = x x {k + 1) + 2s 

= X X ( Xl {k + l),s) , 
ip 2 (Xi (xi(k), s) , u 2 (k, s)) = xi(k) + u 2 (k) + s = x 2 (k + 1) + s 

= X 2 (x 2 (k + l),s) , 
<P3 (X 2 (x 2 (k),s) , ui(k, s)) = (x 2 (k) + s) (wi(fc) + s) 

= x 2 (k)u 1 (k) + s (x 2 (k) + ui(kj) + s 

= x 3 (k + 1) + sxx(k + 1) + 5(s) = X 3 (k + 1) + 6{s) . 

By Theorem 13. II we obtain the following conservation law for the problem: 

2?/;o (xi(k) + x 2 (k)) + 2ipi(k) + tp 2 (k) + tp 3 (k)x 1 (k) = constant. (12) 

Using the information from the discrete time maximum principle, condition 112JI 
is equivalent to 

(xi(k) + x 2 (k)) + 2u 2 (k) — u\{k) — constant . (13) 
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The conservation law (|13fl is a necessary optimality condition. It is trivially 
satisfied choosing the control variables according to: 



Ui(fe) = Xi(fc) , 
U2(k) = --x 2 (k) . 

An extremal is then obtained with the co-state variables given by 

ipl(k) = -u 2 (k) = -x 2 (k) , 
ip 2 (k) = Mi(fc) = xi(k) , 
lfo(*) = 0. 

5 Discrete Calculus of Variations 

We now obtain a discrete Noether's theorem for the problems of the discrete 
time calculus of variations which are quasi-invariant with respect to infinitesimal 
transformations having p parameters, p > 1, up to a difference gauge term. 

5.1 The fundamental Problem 

The fundamental problem in the discrete calculus of variations is a special case of 
our problem (P): r = n; no restrictions on the controls (£1 = R n ); <p{k, x, u) = u. 
The problem is then to determine a finite sequence x{k) £ M. n , k = M, . .. ,M + 
N, x(M) = xm, x(M + N) = xm+n, f° r which the discrete cost function 

M+N-l 

J[x(-)}= L(k,x(k),x(k + 1)) 

k=M 

is extremized. The maximality condition in the Theorem l2.Il implics in this case 
the conditions 

OH 

— (k, x(k),u(k), ifo, Mk + 1)) = , k = M, . . . , M + N - 1 , 
ou 

that is, 

BT 

^{k + 1) = (k, x{k), x(k + 1)) , (14) 

while from the adjoint system one gets 

BT 

yj(k) = ^ —(k,x(k),x(k + l)). (15) 

We note that no abnormal extremals exist for the fundamental problem of the 
discrete calculus of variations: ipo = implies that ip(k + 1) is zero for all 
k = M, . . . , M + N — 1, a possibility excluded by the discrete time maximum 
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principle. So it must be the case that r/'o ^ 0. From l(Ti|) and l|15|). a necessary 
condition to have an extremum is that x(k), k = M, . . . , M + N — 2, must satisfy 
the second-order difference equation 

8L 8L 

— (k + l,x(k + l),x(k + 2)) + — (k,x(k),x(k + l)) = 0. (16) 

Equations l|16|) share resemblances with the continuous Euler-Lagrange equa- 
tions, and are called the discrete Euler-Lagrange equations. 

Definition 5.1. The discrete Lagrangian L (k, x(k), x(k + 1)) is said to be 
quasi-invariant with respect to the infinitesimal p-parameter transformation 
X(k, x,u,s), s — (s%, . . . , s p ), \\s\\ < e, X{k, x, u, 0) = x for all k = M, . . . , M + 
N — 1, x, u € R", up to the difference gauge term $ (k, x(k), x(k + 1), s), if for 
each k = M, . . . , M + N - 2 

L (k, x(k), x(k + 1)) + A$ (k, x(k), x(k + 1), s) + S (k, x(k), x(k + 1), s) 

= L(k,X(k,x(k),x(k + l),s),X(k + l,x(k + l),x(k + 2),s)) , (17) 

where £(-, •, •, •) is a function satisfying Q. 

Corollary 5.1. If L (k,x(k),x(k + 1)) is quasi-invariant with respect to the p- 
parameter transformation X up to the difference gauge term <f>, in the sense of 
Definition \5.1\ then all solutions x(k), k — M, . . . , M + N — 2, of the discrete 
Euler-Lagrange difference equation ()16|l satisfy 



BT B 

— (k- l,x(k - l),x(k)) ■ —X {k,x{k),x{k + l),s) 
ou Bsi 



8 

—$(k,x(k),x(k + l),s) 
8s, 



s=0 

= constant, 



s=0 



5.2 Higher Order Discrete Problems 

Let us now consider the problem of optimizing 

(k,x(k),x(k + 1), . . . ,x(k + m)) , (18) 

k 

where L(k ) is continuously differentiable with respect to all vari- 

ables. This problem is analogous to the continuous problems of the calculus of 
variations for which the Lagrangian L depends on higher-order derivatives. It 
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is easily written in the optimal control form (P). Introducing the notation 

x°(k) = x(k) , 
x 1 {k) = x(k + l), 

x™- 1 ^) = x(k + m-l) , 



one gets: 



u(k) = x(k + m) 



L (k, x°(k), x" 1 " 1 ^), u(k)) — > extr . 

k 

f x°(k + l) = x 1 {k), 
x 1 (k + 1) = x 2 (k), 



x m - 2 (k + l) = x m ~ 1 {k) , 
[x^ik + l) =u(k). 



The Hamiltonian is given by 



' m-2 



H = ^ L(k,x°,...,x m ~ 1 ,u) + Y^ -x^ 1 



3=0 



From the maximality condition 



8T 

^-\k + 1) = (k,x°(k), . . . , x m -\k),u{k)) 



while from the adjoint system 

dL 



^ (k)=^Og^(k,x (k),...,x m - 1 (k),u(k)) , 



dL 



^(k) = Vo^j (k, x°(k), . . .,x m -\k),u(k)) + ^-\k + 1) , 



(19) 

(20) 
(21) 



j = 1, . . . , m — 1. From l|19|) . (|20|l . and (|21|) . we conclude that: similarly to the 
fundamental problem of the calculus of variations, no abnormal extremals exist 
in the higher order case; the equation 

dL 



^2~Q-f (k + m-j,x°(k + m-j),. 



3=0 



V m 1 (k + m-j),u(k + m-j))=0 

(22) 

holds. Going back to the initial notation, (|22H is nothing more than the discrete 
Eulcr-Poisson equation of order 2m for the m-th order discrete problem of the 
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calculus of variations H18|) : 
m dL 

— — r (k + to — j, x(k + m — j), . . . , x(k + 2m — 1 — j),x(k + 2m — j)) = . 



(23) 

Definition 5.2. The discrete Lagrangian L (k, x(k), . . . , x(k + m)) is said to 
be quasi-invariant with respect to the infinitesimal /^-parameter transformation 
X(k,x°,...,x m ,s), s = ( Sl ,...,s p ), \\s\\ < e, X(k,x a ,...,x m ,0) = x" for all k, 
and x J , j — 0, . . . , to, up to the difference gauge term $ (k, x(k), . . . , x{k + to), s), 
if for each k 

L (k, x(k), . . . , x{k + m))+A$ (k, x(k), . . . , x(k + to), s)+5 (k, x(k), . . . , x(k -f to), s) 

= L(k,X (k, x(k), . . . , ar(fc + m), s) , . . . , X (k + m, x(k + to), . . . , x(k + 2m), s)) , 

(24) 

where ^ = 0, i = 1, . . . , p. 

Corollary 5.2. If L (k, x{k), . . . , x{k + m)) is quasi-invariant with respect to 
the p-parameter transformation X up to the difference gauge term <!>, in the 
sense of Definition 15. SI then all solutions x(k) of the discrete Euler-P oisson 
difference equation l|23|) satisfy 



d_ 

dsi 



$ (k, x(k), . . . , x(k + to), s) 



s=0 



m— 1 j 



<9L 



Z Z T~r (/c + J ' ~ ^' T(fc + 3 ~ ■ ■ ■ + z + m)) 



9a; 

3=0 i=0 

• -^~X (k + j, x(k + j), . . . , x(k +j + to), s) 



= constant, (25) 



1, 



In the case to = 1 the discrete Euler-Poisson equation i|23|) reduces to the 
discrete Euler-Lagrange equation (|16|l , and the conservation law (|25|l reduces to 



d 

— $(k,x(k),x(k + l),s) 

OSi 



s=0 



dL d 
+ (*, &(*). a?(fc + 1)) • -^-X (k, x(k), x(k + 1), s) 

or, which is the same, 



constant , 



s=0 



9 

—$(k,x(k),x(k+l),s) 

OSi 



s=0 



r\-r Q 

- — (fc - 1, ar(A - 1), a?(fc)) • ^—X (k, x(k), x(k + 1), s) 

OX L OSi 

This is precisely the conservation law given by Corollary 15. II 



constant . 



s=0 
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